Mechanism of unidirectional emission of ultrahigh Q Whispering Gallery mode in microcavities 
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The mechanism of unidirectional emission of high Q Whispering Gallery mode in deformed circular micro- 
cavities is studied and firstly presented in this paper. In phase space, light in the chaotic sea leaks out the 
cavity through the refraction regions, whose positions are controlled by stable islands. Moreover, the positions 
of fixed points according to the critical line in unstable manifolds mainly determines the light leak out from 
which refraction region and the direction of the emission. By controlling the cavity shape, we can tune the leaky 
regions, as well as the positions of fixed points, to achieve unidirectional emission high Q cavities with narrow 
angular divergence both in high and low refractive index materials. Especially for high index material, almost 
all Gibbous-shaped cavitiess have unidirectional emission. 
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In recent years, being the high potential elements for in- 
tegrated devices, dielectric microcavities have attracted more 
and more attentions. Especial attentions are focused on the 
Whispering Gallery (WG) micro-resonators, such as micro- 
sphere, micro-disk, micro-toroid, etc, where energy can be 
well confined in these rotational symmetrical structures by to- 
tal internal reflections. The unique high quality (Q) factor and 
low mode volume features of WG Modes lead to wild applica- 
tions, ranging from nonlinear optics, low threshold lasers, sen- 
sitive sensors, to cavity quantum electrodynamics Ij]]. How- 
ever, the isotropic emission property of WGMs makes very 
low efficient coupling from free space. Also, the near field 
coupling component is a big obstacle for practical applica- 
tions. To solve the coupling problem, researchers found that 
the microcavity with slightly deformed circular boundary, also 
called asymmetric resonant cavity (ARC), can lead directional 
emission. The ARCs, such as the extensively studied quadru- 
ple lUSIlll Hi and stadium microcavities 1011 H], can be 
efficiently pumping through free space focused beam, which 
makes them have been well used in the realization of strong 
coupling ifioll and laser ifTlll . Besides the easy free space cou- 
pling in the optical experiments, the deformed microcavities 
can also serve a s op en billiards for experimental research on 
quantum chaos lll2lll3ll . 

In order to get more efficient free space excitation and col- 
lection in practical application, people are pursuing single di- 
rectional emission with narrow angular divergence, such as 
the spiral shaped ifTsi [Tsi [T6tl and the rounded isosceles tri- 
angle shaped ITtIi microcavities. However, their low Q fac- 
tors highly limit the applications in low threshold laser and 
cavity quantum electronic dynamics. Although there are sev- 
eral theoretical approaches ifTsl the trade-off between 
the high Q and unidirectional emission still blocks its fur- 
ther development. It is till recently that fabrication-friendly 
limagon-shaped microcavity with both unidirectional emis- 
sion and high Q was theoretically proposed |20] and exper- 
imentally realized by several groups Ii21i, |22|, i23l1 . However, 
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the mechanism to design such cavity is not clear yet. More- 
over, this limagon-shaped microcavity is based on the high 
refractive index material. There is no report for high Q and 
unidirectional emission cavity with low index material. 

In this paper, we will show the mechanism to get unidi- 
rectional emission high Q mode both in the high and low re- 
fractive index ARC. At the beginning, we summarize the nec- 
essary conditions to achieve the high Q unidirectional emis- 
sion ARC. First, the cavity boundary should be continuous, 
smooth, and slightly deformed from a circle to support the 
WG-like modes. In such deformed microcavity, the wave is 
mainly localized high above the critical line and well confined 
in the cavity to support high Q. The light may leak out when 
the ray is lower than critical line in phase space. In this case, 
the rays follow the unstable manifolds, cross the critical line, 
and perform the directional emission tangential to the cavity 
boundary. 

Second, the boundary shape should have less symme- 
try. For planar deformed cavity with two or more axial 
symmetries, there are at least two pairs of symmetric tan- 
gential emission spots on the boundary because they sup- 
port both clockwise and counterclockwise traveling modes 
I|2i,i3|,i4i,^^ 7,_8,^. The only way to achieve unidirectional 
emission is to design the cavity with at most one symmetric 
axis and make the emission direction along the axis direction. 

In the cavity with continuous, smooth, and slightly de- 
formed circular boundary, there are stable and unstable or- 
bits, corresponding to the stable islands and fixed points of 
unstable manifolds in phase space. The mechanism for the 
unidirectional emission in the single-axis cavity is: The light 
only along the manifolds may leak out the cavity. Isolated 
by the stable islands, the manifolds have several leaky regions 
when they cross the critical refractive line. By controlling the 
cavity shape, we can tune the leaky regions, as well as the po- 
sition of fixed points according to the critical line. The light 
through the lower fixed point has higher probability to leak out 
than through the higher point. The parallel leaky lights along 
the axis direction shows unidirectional emission from the high 
Q ARC. Based on this mechanism, we succeed in designing 
the unidirectional high Q ARC both in high and low refrac- 
tive index materials. We find that almost all Gibbous-shape 
ARCs have unidirectional emission for high refractive index 
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FIG. 1 : (color online) (a) The real space illustration of the rays reflec- 
tion inside HQHC deformed microcavity with e = 0.11. The dashed 
line outside is the circle boundary. The blue lines and red dash lines 
are the stable and unstable period-3 orbits, (b) The calculated field 
spatial distribution of TM polarized fundamental WG-like mode in 
HQHC, with the wavenumber in vacuum kr ~ 26.15 and quality 
factor Q = 3.3 x 10"". (c) The phase space structure of HQHC. The 
blue circles and red diamonds are the stable islands and fixed points, 
corresponding to the stable and unstable period-3 orbits in (a). The 
red line is the critical refraction line with sin x = 1/?t-, and the green 
points are the unstable manifolds from the unstable period-3 orbits. 
Magenta points is the collection records of leaky rays location from 
the original rays upon sin a:: > 0.7. It shows that most rays leak out 
from this refraction region, (d) the far field patterns for TM polar- 
ized ray (the red line, lifted by 0.4) and wave simulation (blue line 
and green line for even and odd parity mode respectively). 

material. By precisely controlling the boundary shape, we can 
achieve far field emission with very narrow divergence angle. 

At first, we consider a boundary shape of Half-Quadruple- 
Half-Circle (HQHC) to illustrate the unidirectional emission 
in high refractive index cavity, which have been demonstrated 
in silica microsphere |24j]. The slightly deformed circular 
boundary shape, as shown in Fig. 1(a), can be read as 

, . _ ( i?o, cos(/> < 0, 

" \ i?o (1 - ecos20) , cos > 0, 

where (p is the polar angle according to the symmetric X-axis 
and e is the deformation factor. In the limit of X, the light 
in cavity can be semiclassically treated as ray. Although wave 
nature of light is important in microcavity, the ray dynamics 



is still an efficient tool to understand and analyze the emission 
properties of ARCs 1,22), ,23t, |25|, i26t (27^ |28t, |29(, i3ai . 

The sequence of ray reflection inside the cavity in real 
space can be represented in the phase space as Poincare Sur- 
face of Section (SOS). Each ray reflection on the boundary is 
recorded in Birkhoff coordinates ll3Tll as (</>, sinx), where (j) 
and X denote the polar angle of reflection position on bound- 
ary and the incidence angle, respectively. By setting rays ini- 
tially random on boundary (here we only consider the coun- 
terclockwise propagation rays, the clockwise SOS can be ob- 
tained by symmetric transformation), we can get the SOS by 
treating the cavity as billiard lfl2ll . which is shown in Fig. 1(c). 

In SOS, the phase space structure can be divided into two 
parts, the chaotic sea and the stable islands. The two regions 
are isolated to each other The ray in chaotic sea [black points 
in Fig. 1(c)] has chaotic motions in real space, while the ray in 
islands (in blue) is conserved and will never move into chaotic 
sea. In the real space, the stable islands are corresponding 
to the stable (blue solid lines) period-3 orbits 1.31.1 , as shown 
in Fig. 1(a). The hght in the chaotic sea will run along the 
unstable manifolds [26] [green curve in Fig. 1(c)] and leak 
out when it crosses critical line, which is denoted by the red 
line in Fig. 1(c). The critical line shows the critical incidence 
angle for total internal reflection sinx = i/n. We set the 
refractive index n — 3.3 in our case. 

In this slightly deformed circular cavity, the light is mainly 
localized upon the critical line in phase space. In the real 
space, the light runs in the unstable periodic orbits or high 
Q WG-Iike mode, as shown in Fig. 1(b) SIMIll- Initially, 
we set rays uniformly distributed in the region sinx > 0.7, 
which can be considered as energy of high Q modes dynamics 
localized above the critical line. Rays will be reflected on the 
boundary in sequence until they hit the refraction regions. In 
the SOS, the refraction regions are the parts of chaotic sea be- 
low the critical line (sinx < ^/n). As the manifolds enclose 
the stable islands, the locations of the stable islands determine 
the positions of these refraction regions. In our case, the cen- 
ter two stable islands are lower than the third island. So, there 
are two refraction regions near i/i w 7r/2 and 37r/2 at bound- 
ary, corresponding to two directional emissions with far field 
angles of 180 and 360 degrees. 

In order to get the unidirectional emission, we need to ana- 
lyze the relative positions of the three fixed points in the man- 
ifolds, as shown with red diamonds in Fig. 1(c). It is corre- 
sponding to the unstable (red dashed lines) period-3 orbit lIsTll 
in Fig. 1(a). They are sandwiched by the stable islands. All 
lights following the manifolds will transit through/near these 
three points. However, the fixed point (marked as 1) before re- 
fraction region near (j) w 7r/2 is lower than the point (marked 
as 2) before the other refraction region near 4> ~ 37r/2. Also, 
it is closer to the refraction region. So the manifolds here 
are shorter and steeper to cross the critical line. As a result, 
the light through/near fixed point 1 has the higher probabil- 
ity to enter the fraction region and leak out than the other 
one. It will give rise to the most energy leaky here and the 
directional emission with far field angle of 180 degrees. From 
the symmetry, the clockwise propagation rays leak out around 
(f) w 37r/2 and show the same directional far field emission. 
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Therefore, the leaky beams from the counterclockwise and 
clockwise propagations are parallel and finally show unidi- 
rectional emission from this HQHC cavity. 

With the numerical solution of Maxwell equations through 
boundary elements method [32., ,33]. we can observe the high 
Q WG-like mode in this HQHC microcavity. Fig. 1(b) shows 
the electric field intensity distribution of WG-like transverse 
magnetic (TM) polarized mode with radiation quantum num- 
ber q — 1 (the transverse electric polarized modes have the 
similar properties). The near field pattern shown tangent light 
emission from « 7r/2 and 3tt/2, corresponding to the two 
refraction regions. Obviously directional emission beams for 
(counter)clockwise propagation are along the X-axis. 

The far field distribution of the Odd and Even polarity WG- 
mode is shown in Fig. 2(d), as well as the ray-wave corre- 
spondence with ray dynamics. In the ray simulation, with the 
weighted refractive coefficient from Fresnel law, we can get 
the far field amplitude by summing up the all refractions, as 
shown in Fig. 2(d) with the red curve. In order to get clear 
illustration, the curve is lifted by 0.4. The far field pattern 
shows the unidirectional emission at 180 degree with diver- 
gence about 30 degrees. 

As discussed above, the refraction regions and fixed points 
in unstable manifolds play the most important roles to deter- 
mine the cavity directional emission. So, we can tune the 
cavity shape to adjust the locations of stable islands and the 
positions of fixed points to control the refraction regions and 
the emission direction, simultaneously. To achieve the unidi- 
rectional emission of high Q modes, we can design the X-axis 
symmetric cavity and lead the unidirectional emission along 
the axis. In general, the boundary of X-axis symmetric shape 
can be expressed as 



Ro J2 o.i cos* 

i?0 cos' 



COS 4> > 
COS (h <0 



(2) 



By setting ao = &o = 1, and ai = 6i = 0, the norm direction 
is continuous and the cavity boundary is smooth. In addition, 
we should keep the boundary slow varying, so simply we cut 
off the high order terms, only keep 02 (&2) and 03(63) nonzero 
lf34ll . Also, to break Y-axis symmetry, it needs a2 ^ 62 or/and 
0-3 7^ ^^3- Moreover, we can set 62 = 63 = and 02 + 03 < 
to form Gibbous shape for further simplification. 

For the this kind of Gibbous shape cavity with high re- 
fractive index material, there will always exist stable and un- 
stable period-3 orbits with vertexes near (0, 2tt/3, 47r/3) and 
(tt/S, tt, 57r/3), respectively. By randomly set the parameters 
fli, we can always get similar phase space structures to Fig. 
1(c) in any Gibbous-shaped cavity. That is because, without 
Y-axis symmetric, the two sets of period-3 orbits in Gibbous 
shape cavities always have smaller convex angles near 7r/2 
in real space. Correspondingly, the fixed points of manifolds 
and stable islands are lower in phase space, which make the 
refraction regions fixed at about (j) k, n/2 and 37r/2 and uni- 
directional emission from the first refraction region. We car- 
ried out ray and wave simulations on lots of Gibbous shape 
cavities. Results indicate the single emission with far field di- 
vergence angle ranges from 20 to 50 degrees, confirming to 
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FIG. 2: (color online)(a) (b) the SOS of two cavities (a) HQHC 
e = 0.05 and (b) the cavity shape support stable rectangle orbit and 
unidirectional emission designed in the text, (c) the rectangle and di- 
amond shape period-4 orbits in the X-symmetric cavity, (d) The near 
field distribution of the unidirectional emission with kr ~ 52.5284. 
(e) The far field pattern of the cavity, red line is the ray simulation 
result (lifted by 1), and the blue line is corresponding to the result in 
(d). 



our conjecture from the orbits. We find the cavity with 02 ~ 
0.0486, 03 — —0.1258 gives good unidirectional far field pat- 
tern. The wave simulation gives a high Q{Q ~ 5.78 x 10^) 
TM polarized WG mode with kr « 25.1842, and the far field 
divergence is only about 24 degrees, which is much smaller, 
comparing to 40 degrees reported in lima§on cavity |[20ll . 

Thus, we have shown how to get unidirectional emission 
in the Gibbous-shaped cavity. The limagon cavity ll20ll . as 
well as the cavity in Song's experiments [211], is a generalized 
Gibbous-shaped cavity. Their mechanisms for unidirectional 
emission can also be explained with our approach. 

Now, we will generalize the above mechanism for unidirec- 
tional emission to low refractive index microcavities. Here we 
take the silica {n — 1.45) for example. Similar to the period- 
3 orbits in high index cavities, the silica cavity's directional 
emission is basically influenced by period-4 orbits. The for- 
mer experimental study on HQHC micro-sphere showed that 
the emission direction could be reduced to 2 nearly perpen- 
dicular direction [24]. The SOS of such HQHC is present in 
Fig. 2(a), where the stable period-4 islands [the red diamond 
orbit in Fig. 2(c)] around the critical line prevent the rays uni- 
directional emission around </) = 7r/2. 

A strategy to solve the problem in this low index cavity is 
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to find a shape with stable islands away from 7r/2. That is to 
make the green rectangle period-4 orbit be stable and the red 
diamond orbit be unstable. Corresponding SOS is presented 
in Fig. 2(b). Rays can go cross the critical line near 7r/2. 
Similar to the analysis in high index cavity, the conditions for 
unidirectional emission in silica microcavity is: (i) The cavity 
supports the stable rectangle orbit, (ii) The stable island near 
(f) — 3tt/A and 57r/4 is lower, and the fixed point (the red 
triangle in SOS) is lower near (f> — tt/2. 

In the numerical simulation of ray dynamics, we randomly 
choose parameters to set the boundary shape, judge that the 
rectangle orbit stable or not, and compare the islands posi- 
tion in SOS to fulfill the condition (ii). We can easily find a 
good cavity shape with 02 = —0.1329,03 — 0.0948,52 — 
—0.0642, 63 = —0.0224 satisfies the above conditions, and 
gives good unidirectional far field pattern. Fig. 2(b) and 
Fig. 2(d) shows the SOS and field distribution of TM mode. 
Fig. 2(e) is the corresponding far field pattern compares to the 
ray simulation. As we expect, the Q factor is up to 2.75 x 10^ 
and divergence is only about 30 degree. 

Thanks to the ray-wave correspondence, an efficient route 
to design a unidirectional emission cavity shape could be: us- 
ing the ray dynamics to give the SOS of the selective boundary 
shape, predicting the emission properties, and carrying out the 
wave simulation to confirm it. Here, we only gives examples 
on the typically materials, such as semiconductor and silicon 
(n « 3.3) and silica glass {n « 1.45). Other popular solid 
materials for photonics have similar refractive index. There 
are also many materials not at this range, such as some doped 
glass with the refractive index around 2.0. We have also suc- 
cessfully designed the cavity shapes to achieve high Q and 
unidirectional emission (not presented here). By adjusting the 
periodic-3 orbits or periodic -4 orbits, we can get appropriate 
emission direction and expect the unidirectional emission in 
the materials with refractive index range from 1.4 to 4. 

There are also some deviations between the simulations 
with ray and wave, especially when kr is small. The actions 



of light could not totally be described by rays. Better ray- 
wave correspondence could be expected in larger cavity ll23ll . 
In smaller cavity, wave properties should be included, such as, 
the diffractions and interference. In microcavities, researches 
have illustrated that some modifications should be included 
in the ray dvnamics ||35n. such as the corrected Fresnel's law 
at curve interface ll36ll . the Goos-Hanchen shift ifssll and the 
Fresnel Filtering effect iIstIi . 

In conclusion, we examined and presented the mechanism 
for high Q unidirectional emission WG modes in micro- 
cavities. Based on the necessary condition for the continu- 
ous and single axis symmetric boundary shape, we can well 
control the directional emission from ARC by setting appro- 
priate boundary shape to adjust the stable islands and fixed 
points of unstable orbits in phase space. With the assistance 
of ray dynamics, it is easy to design fabrication-friendly sim- 
ple cavity boundary shape to achieve high Q and unidirec- 
tional emission in different materials with the refractive index 
range from 1.4 to 4. We expect our approaches presented here 
to discuss the light in cavities can be applied to the study of 
the chaotic transport in two dimensional phase space |39], as 
well as other opening non-integrable systems, such as quan- 
tum dots and nano-structures. 
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